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Abstract 

We investigate the space-time evolution of ultrarelativistic Au-Au collisions at full RHIC energy 
using a schematic model of the expansion. Assuming a thermally equilibrated system, we can 
adjust the essential scale parameters of this model such that the measured transverse momentum 
spectra and Hanbury-Brown Twiss (HBT) correlation parameters are well described. We find that 
the experimental data strongly constrain the dynamics of the evolution of the emission source 
although hadronic observables for the most part reflect the final breakup of the system. 
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I. INTRODUCTION 



Numerical simulations of finite temperature Quantum Chromodynamics (QCD) on discrete 
space-time lattices suggest that the theory exhibits a transition to a new state of matter, the 
quark-gluon plasma (QGP), at temperatures of the order of 170 MeV Experimentally, 
hadronic matter at this temperature can only be created for a short time in ultrarealtivistic 
heavy-ion collisions, and currently there are ongoing efforts from both experiment and theory 
to find conclusive evidence for the creation of the QGP and ultimately to study its properties. 
A wealth of data on hadronic single particle distributions and two particle correlations has 
been assembled so far. 

In the interpretation of experimental data of ultrarelativistic heavy-ion collisions, however, 
one is often faced with the challenge to disentangle signals of new physics indicating the 
production of a QGP from known hadronic effects. This question can only be reliably 
addressed if the space-time evolution of the fireball created in the collision is sufficiently 
known. In this note, we present an attempt to determine the evolution of the bulk hadronic 
matter by fitting the essential scales of a schematic model for the fireball expansion dynamics 
to the measured hadronic data. The resulting scenario can then be used in the calculation 
of other observables (not used in the fit) which are sensitive to the bulk matter expansion, 
such as the emission of dileptons and photons or jet quenching, thus reducing or eliminating 
the inherent ambiguity in the interpretation mentioned above. 



II. THE MODEL 



The main assumption for the model is that an equilibrated system is formed a short time 
tq after the onset of the collision. Furthermore, we assume that this thermal fireball subse- 
quently expands isentropically until the mean free path of particles exceeds (at a timescale 
Tf) the dimensions of the system and particles move without significant interaction to the 
detector. In addition to this final breakup (freeze-out) of the fireball, particles are emit- 
ted throughout the expansion period whenever they cross the boundary of the thermalized 
fireball matter. 

For simplicity we restrict the discussion to a system exhibiting radial symmetry around the 
beam (z)-axis corresponding to a central {b = 0) collision. For the entropy density at a given 
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proper time we make the ansatz 



s(r, rjs, r) = NR{r, r) ■ H{t]s, r) 



(1) 



with T the proper time as measured in a frame co-moving with a given volume element 



shape of the distribution and a normalization factor. We use Woods-Saxon distributions 



to describe the shapes for a given r. Thus, the ingredients of the model are the skin 
thickness parameters d^s and ?7„s and the parametrizations of the expansion of the spatial 
extensions Rc{t),Hc{t) as a function of proper time. For simplicity, we assume for the 
moment a radially non-relativistic expansion and constant acceleration, therefore we find 
Rc{t) = Rq + ^T^- Hc{t) is obtained by integrating forward in r a trajectory originating 
from the collision center which is characterized by a rapidity i]c{t) = i]q + a^^r with rj^ = 
atanh f ^ the longitudinal expansion velocity for that trajectory. Since the relation between 
proper time as measured in the co-moving frame and lab time is determined by the rapidity 
at a given time, the resulting integral is in general non-trivial and solved numerically (see [2] 
for details). Rq is determined in overlap calculations using Glauber theory, the initial size of 
the rapidity interval occupied by the fireball matter, rjo is a free parameter and we choose to 
use the transverse velocity v{_ = a±Tf and rapidity at decoupling proper time f]^ = f]o + a^Tf 
as parameters. Thus, specifying riQ,rif,v{^ and Tf sets the essential scales of the spacetime 
evolution and rf^s and r^ws specify the detailed distribution of entropy density. For simplicity, 
we do not discuss a (possible) time dependence of the shapes (e.g. parametrized by d^siT) 
and ?7„s(t)) at this point. 

We require that the parameter set {rjf, ?7„s) reproduces the experimentally observed rapidity 
distribution of particles but we do not make any assumptions about the initial rapidity 
interval characterized by rjo. This allows for the possibility of accelerated longitudinal ex- 
pansion and implies in general r] ^ r]s- Here, = | In ^^3^ denotes the longitudinal rapidity 
of a volume element moving with momentum p'^. We require that the longitudinal flow pro- 
file is such that an initially homogeneous distribution remains homogeneous. For a general 
accelerated motion, there is no simple analjd;ical expression for the spacetime position of 



and Tjs = ^ ln( 



-) the spacetime rapidity and R{r, r), H{ris, r) two functions describing the 
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sheets of given r or the flow profile corresponding to r = y/t"^ — z"^ and r] = rjg valid for 



the non-accelerated case. In P] however, we have investigated such an expansion pattern 
with a constant acceleration and argued that for rapidities < 4 the acceleration leads to an 
approximately linear relation rj ^ (rjg at given r and that sheets of constant proper time 
are approximately hyperbolae as in the non-accelerated case. This mismatch between rj and 
7]s leads to additional Lorentz contraction factors in volume integrals at given r, hence the 
longitudinal extension of matter on a sheet of given r in the interval — 7^^''°°* < rjs < rif""^^ 
must then be calculated as 

Hr) « 2r-'°" (3) 

In the following, we use these results in our computations whenever rjf 7^ rjQ. For rjf = rjo, 
the model reduces to the well-known expressions of the Bjorken expansion scenario, e.g. 



L(r) = 2riQT with r = vt^ 



z 
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For transverse flow we assume a linear relation between radius r and transverse rapidity 
p = atanh v±{t) = r / Rc{t) ■ Pc{t) with Pc{t) = atanh a±T. For the net baryon density 
inside the fireball matter we assume a transverse distribution (apart from a normalization 
factor) given by -R(?vr), but its longitudinal distribution we parametrize such as to describe 
the measured data 

We proceed by specifying the equation of state (EoS) of the thermalized matter. In the QGP 
phase, we use an equation of state based on a quasiparticle interpretation of lattice QCD 
data (see 0]). In the hadronic phase, we adopt the picture of subsequent chemical freeze-out 
(at the transition temperature Tc) and thermal freeze-out (at breakup temperature Tp) and 
consequently use a resonance gas EoS which depends on the local net baryon density. The 
reason for this is that a finite baryochemical potential pb leads to an increased number 
of heavy resonances at the phase transition point, and decay pions from these resonance 
decays lead in turn to a finite pion chemical potential p.,^ in the late evolution phases which 
implies an overpopulation of pion phase space and faster cooling as compared to a scenario 
in chemical equilibrium. We calculate the pion chemical potential as a function of the local 
baryon and entropy density using the statistical hadronization framework outlined in [6] . We 
find Pt^ to be small (0(30 MeV)) in the midrapidity region at RHIC, however the corrections 
to the chemically equilibrated case are important. 

With the help of the EoS, we can find the local temperature T{rjs, r, r) of a volume element 



from its entropy density s{r]s, r, r) and net baryon density psiVs, t)- We calculate particle 
emission throughout the whole lifetime of the fireball by selecting a freeze-out temperature 
Tf, finding the hypersurface characterized by T{ris,r,T) = Tj and evaluating the Cooper- 
Frye formula 



E- 
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da^p^ exp 



Tf 



d'^xS{x,p) 



(4) 



(27r)3 

with the momentum of the emitted particle and g its degeneracy factor. Note that the 
factor da^j^p^ contains the spacetime rapidity rjg and the factor p^u^ the rapidity rj. Since 
these are in general not the same in our model, the analytic expressions valid for a boost- 
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FIG. 1: Left panel: Time evolution of the average fireball temperature for the fitted set of pa- 
rameters (solid line) and a model assuming boost-invariant longitudinal expansion (dotted). Right 
panel: Measured transverse momentum spectra for vr^ (circles), (squares) and p (diamonds) 
as compared to the model results. 

In our case, the parameter tj denotes the freeze-out time for the last volume element to 
reach the temperature Tj. Volume elements freeze out throughout the whole evolution time 
(and are boosted with a time- and position dependent flow velocity) whenever they cross 
the Cooper-Frye hypersurface. Therefore the transverse expansion parameter v^^ also does 
not necessarily reflect the typical transverse velocity of a volume element since freeze-out 
may occur earlier than Tf and not at the radius r = R(t). 

nn 

Using this emission function, we calculate the HBT parameters as [3, 



Z\2\ 



'z2\ 



(5) 
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with x^ = Xfj, — {Xfj,) and 



Jd^xSix,K) 



In order to parametrize the spacetime evolution of a central 200 AGeV Au-Au collision 
at RHIC we fit the remaining set of parameters dwH, Vn, Tf and f;[ to the experimentally 
obtained single particle transverse momentum spectra jlo| and HBT parameters For 
the description of the HBT correlation parameters which are measured for 30% central 
collisions, we scale down the entropy content and initial overlap radius guided by overlap 
calculations and neglect the angular asymmetry. 



III. RESULTS AND DISCUSSION 

We find that the choice of parameters Tf = 110 MeV, d^s < 1-0 fm, 7]o = 1.8 and = 0.67 is 
able to give a good description of the data. The resulting average cooling curve (computed by 
averaging the entropy density over the fireball volume at a given proper time and determining 
the corresponding temperature, where we define the fireball volume at given r as the 3- 
volume bounded by the Cooper-Frye surface) and the transverse momentum spectra for 
7i~, K~ and nucleons are shown in Fig. the HBT correlation parameters in Fig. |21 
The fit apparently misses the low pt part of the transverse momentum spectra for the heavier 
particles but describes the 7r+7r+ correlation radii well, with the exception of the low pt part 
of i?out where the calculation lies somewhat above the data. 

The rather steep falloff of Rsidei -Rout as a function of kt favours a large amount of fiow. How- 
ever, demanding simultaneous agreement with the slope of the momentum spectra implies 
that large transverse fiow has to be accompanied by low freeze-out temperatures. There- 
fore, the volume at freeze-out has to be large in order to reach small entropy densities. The 
inclusion of a finite pion chemical potential (originating from resonance decays in chemical 
non-equilibrium) is crucial to reduce the volume corresponding to a given average temper- 
ature. Since the radial expansion is rather constrained by the normalization of the HBT 
parameters, a large volume implies sizeable longitudinal extension and hence a long-lived 
system. Such a long lifetime in combination with a standard boost-invariant longitudinal 
expansion leads to large values of -Riong which are clearly incompatible with the data. A size- 
able longitudinal compression and re-expansion however leads to a reasonable description of 
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FIG. 2: The HBT correlation parameters R^e, Ro^t, Riong and the ratio i?out/^side in the model 



calculation as compared to PHENIX data 



-Riong as well. 

It is not possible to identify a single feature of the model as being responsible for the good 
quality of the fit but one can illustrate certain trends. Neglecting the additional cooling 
caused by a finite value of fi-,,, it is still possible to describe the spectra and the slope of the 
HBT data well for a low freeze-out temperature, however a larger volume is needed to get 
to this temperature and therefore the normalization of the correlation radii is systematically 
above the data (see Fig. El left panel). 

On the other hand, retaining the idea of a boost-invariant scenario without longitudinal 
acceleration implies a short lifetime of the system as apparent from the estimate -Riong = 
Tf{Tf /rritY/'^ based on such an expansion pattern. Such a short lifetime can be achieved 
by a large value of Tj. The transverse expansion velocity can then be adjusted to fit the 
mt spectra, but the resulting solution is characterized by small fiow velocities which leave 
both normalization and falloff of -Rgide in disagreement with the data (see Fig. El left panel). 
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FIG. 3: Left panel: -Rside in a scenario with vanishing pion chemical potential [i-,^ (dotted line) and a 
difference scenario assuming Ty = 140 MeV (dashed line) as compared to the standard calculation 
(solid) Right panel: iiiong in a scenario with = 110 MeV but assuming no longitudinal acceler- 
ation (dotted) compared with the standard scenario (solid). Note that the transverse momentum 
spectra are well described by all of the scenarios shown in this figure. 

Retaining a low freeze-out temperature in a boost-invariant scenario allows to describe the 
transverse spectra and correlations well but significantly overpredicts -Riong (see Fig. El right 
panel) . 

Concluding, we find that the mt spectra and HBT parameters measured at RHIC can be 
simultaneously described assuming a scenario with small freeze-out temperature Tp ^ 110 
MeV and a sizeable initial longitudinal compression and re-expansion. The data strongly 
constrain alternative scenarios. An upcoming publication will investigate the dependence 
of the results on the different model parameters in more detail. The calculation of further 
observables reflecting different properties of the fireball expansion such as the emission of 
electromagnetic probes or the suppression of high "pt jets within the same framework will 
help to confirm or disprove the outlined scenario. 
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